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ABSTRACT 

We explore the predictions of the standard hierarchical clustering scenario of 
galaxy formation, regarding the numbers and metallicities of PopIII stars likely to 
be found within our Galaxy today. By PopIII we shall be referring to stars formed at 
large redshift (z > 4), with low metallicities ([Z/Zq] < —2.5) and in small systems 
(total mass ^ 2x 10 8 M Q ) that are extremely sensitive to stellar feedback, and which 
through a prescribed merging history (Lacey & Cole 1993) end up becoming part of the 
Milky Way today. An analytic, extended Press-Schechter formalism is used to get the 
mass functions of halos which will host PopIII stars at a given redshift, and which will 
end up in Milky Way sized systems today. Each of these is modeled as a mini galaxy, 
with a detailed treatment of the dark halo structure, angular momentum distribution, 
final gas temperature and disk instabilities, all of which determine the fraction of the 
baryons which are subject to star formation. Use of new primordial metallicity stellar 
evolutionary models allows us to trace the history of the stars formed, give accurate 
estimates of their expected numbers today, and their location in L/ ' L@ vs. T/K HR 
diagrams. A first comparison with observational data suggests that the IMF of the 
first stars was increasingly high mass weighted towards high redshifts, levelling off at 
z 9 at a characteristic stellar mass scale m s = 10 — 15M Q . 

Key words: galaxies: formation - galaxies: evolution - Galaxy: formation - Galaxy: 
abundances - stars: luminosity function, mass function 



1 INTRODUCTION 

The formation and the properties of the first stars in the uni- 
verse are probably among the most fascinating problems in 
present-day cosmology. Numerical simulations and analyti- 
cal studies in a cosmological context (Tegmark et al. 1997, 
Abel et al. 1998, Bromm, Coppi & Larson 2000, Machacek, 
Bryan & Abel 2000) have reached a consensus on the fact 
that the fragments potentially leading to star formation were 
already available inside very small dark matter halos at red- 
shift z w 30 in most cosmological models. After that point 
the requirements concerning the dynamic range and num- 
ber of physical processes to be treated become so demand- 
ing that one has to turn to simplified and/or lower dimen- 
sional studies. Although considerably refined, these studies 
(Omukai & Nishi 1999, Susa & Umemura 2000, Nakamura & 
Umemura 2000, Ripamonti et al. 2000) have not yet reached 
an agreement on the mass of the stellar seed (core) and on 
the accretion rate by which the former grows into a star. 
Determining if and when the accretion is quenched seems to 
hold the key for the understanding of the properties of the 
first stars. However, the large variety of physical effects to be 



considered (for a review see Ferrara 2000) again complicate 
the problem tremendously. 

On general grounds some authors have speculated that 
the IMF of the first stars can be skewed towards high stel- 
lar masses (Larson 1998) or even be bimodal (Nakamura 
& Umemura 2000). One of the consequences of this as- 
sumption is an increase of the relative number of mas- 
sive, short-lived stars. These stars will soon end their life 
as supernovae, leaving either a compact remnant (a neu- 
tron stars or a black hole) or nothing at all (for masses 
above « 12OM0, see Umeda et al. 2000 and Fryer, Woosley 
& Heger 2000), and possibly releasing detectable bursts of 
gravitational waves (Schneider et al. 2000). In turn, this 
implies a strong decrease of the amount of zero- or very 
low-metallicity stars presumably included in the present-day 
Milky Way (MW), according to the hierarchical paradigm. 
Given the above discussion it appears that number counts 
of very low-metallicity MW stars can not only provide evi- 
dences for "living fossils", i.e. stars formed immediately af- 
ter the end of Dark Ages, but also probe the primordial IMF 
and hierarchical models to cosmic epochs (z !t 6) currently 
unreachable to essentially all experiments. 
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We propose here that the first stars are born inside the 
so-called PopIII objects (Haiman & Loeb 1996, Gnedin & 
Ostriker 1997, Ciardi et al. 2000) i.e. systems with total 
mass ^ 2 x 10 s M©, which strongly rely on molecular hydro- 
gen cooling in order to collapse, given their low virial tem- 
peratures; in short, we define as PopIII stars those formed 
inside such mini galaxies. As Ciardi et al. have pointed out, 
these objects are extremely fragile to the energy release by 
their own supernovae; as a result, their gas is quickly evac- 
uated and they can only witness a single burst of star for- 
mation. The remnant is a "naked stellar object" (Ciardi et 
al. 2000), i.e. a tiny agglomerate of long lived low-mass stars, 
essentially retaining the metallicity of the gas from which 
they formed. These are the stars that we are considering as 
candidates for the very low-metallicity stars observed in the 
MW halo. 

Our paper is structured as follows: the next section 
presents the Press-Schechter formalism used to calculate the 
mass functions of MW sub-galactic clumps at different red- 
shifts. In Section 3 we describe the treatment of PopIII halos 
as mini galaxies to estimate the number of stars each con- 
tributes to the final MW. Section 4 gives a comparison with 
recent observational results of the numbers and metallicity 
distribution of extremely low metallicity halo stars. Finally, 
a brief summary is given in Section 5. 



2 BUILD UP OF THE MILKY WAY 

We use the extended Press-Schechter formalism to calculate 
the conditional probability that a virialized halo of mass 
Mi at z\ will become part of a more massive halo of mass 
Mo at a later time, zo- This allows us to calculate the mass 
functions of the fragments which through mergers will end 
up as part of a MW sized system today. Following Nusser & 
Sheth (1999), we identify the mass contained in progenitors 
of mass between Mi and (Mi + dM\) with: 



AM = 



M 5 c {zi — zp) 
2^ (Si ~ So)*/- 



: exp 
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In the above expression \/~Si is the rms density fluctuation in 
a top hat window function of radius (SMi/inpo) 1 ^ 3 , and 5 C 
the critical overdensity for collapse, with po = 3£ImH% /8irG 
being the present mean mass density of the universe. Both of 
the above quantities were calculated using a ACDM power 
spectrum normalized to erg = 0.91, in a Qa = 0.7 flat uni- 
verse (f?M = 0.3), with Ho = 65 km s~ Mpc~ , which we 
used throughout, e.g. Percival & Miller(1999). A baryonic 
fraction of 0.05 of the total halo mass was assumed for all 
halos. 

Current direct estimates of the mass of the MW rely 
on distant galactic satellites used as probes of the mass in- 
ternal to their current positions, using their line of sight 
velocities and reasonable assumptions on their orbits. These 
methods give values for the total mass of the MW of around 
1 x 10 12 M Q , e.g. Kulessa & Lynden-Bell (1992), Binney & 
Merrifield (1998). However, those estimates yield only lower 
limits, as the total extent of the Galactic dark halo could 
be greater than the current positions of suitably observable 
satellites. Consistently with the Press-Schechter approach 
we are using, we estimate the total mass of the MW as the 
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Figure 1. Mass functions of halos which at z = end up as 
part of a 2 X 1O 12 M0 system, i.e. similar to the MW, at different 
redshifts. The dotted lines indicate the mass range where PopIII 
stars are expected to form. 



total mass contained within a sphere of radius r2oo, the ra- 
dius inside which the average density is 200 times the present 
background density, P200 = 200 x po- This in fact is the only 
meaningful definition of "galaxy" within the Press-Schechter 
approach. The circular velocity at this radius would be 
around Vmax = v2 x V200, where Vmax = 220 km s" 1 and 
the factor of v2 between V ma x and V200 is deduced from 
CDM simulations e.g. Avila-Reese et al. (1999), Steinmetz 
& Navarro (1999). This gives: 



: 200 



/3X 1 / 2 
Mmw = I -; — 



(2) 



which for the adopted cosmological model gives 2.6 x 
1O 12 M0, somewhat above the observational lower limits. In 
what follows we shall use Mmw = 2.0 x 1O 12 M0, within the 
uncertainties associated with going from Vmax to V200, and 
in determining Vmax- 

Fig. 1 shows the mass of sub-halos which at z — find 
themselves included in MW sized halos, as a function of sub- 
halo mass, at different redshifts; the dotted lines show the 
limits within which halos at the redshifts shown will form 
PopIII halos. 

We associate a maximum metallicity to the stars to be 
formed and contained within these halos through the mean 
mass weighted metallicity of the universe at each redshift, 
as calculated by Gnedin & Ostriker (1997) using fully hy- 
drodynamic simulations. As a result of the merger process, 
these halos will probably contain a few stars with yet lower 
metallicities formed in halos of the preceding hierarchies; in 
this sense we take the metallicities of Gnedin & Ostriker 
(1997) as an upper limit for the stars of these halos. 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological Origin of the Lowest Metallicity Halo Stars 3 



-4.0 



log(Z/Z s ) 

-3.7 -3.46 -3.1 -2.7 -2.4 




Redshift 

Figure 2. Total mass in PopIII halos that ends up in a MW sized 
galaxy today, as a function of redshift (solid line). The dotted 
curve shows the typical mass of these PopIII halos. The upper 
axis shows the average metallicity of the IGM at that redshift. 



The lower mass limit, M cr u, is taken from Tegmark et 
al. (1997) (see however Fuller & Couchman 2000). Systems 
having virial temperatures below this curve cannot cool in a 
Hubble time. The upper limit, Mt y , corresponds to the up- 
per mass at which the stellar feedback due to supernova can 
expel the baryonic component entirely, thus quenching fur- 
ther star formation. Following Ciardi et al. (2000), we define 
this as a "blowaway", and we adopt their determination for 
the evolution of Mb y . Thus, larger systems where self enrich- 
ment can continue and stars of metallicities higher than the 
Gnedin & Ostriker (1997) limit will be formed are excluded. 
However, given the sharply decreasing mass functions, our 
results are rather insensitive to this upper limit. These two 
limits change slightly with redshift, the ones shown are typ- 
ical for the redshifts displayed, although the more detailed 
redshift dependent estimates discussed above were used. 

Fig. 2 shows the integral of the PopIII halos mass func- 
tions between the selected limits, M tota i, as a function of 
redshift. The upper axis shows the mean mass weighted 
metallicity corresponding to each redshift, discussed above. 
Fig. 2 also gives the typical mass of a single PopIII halo, 
Mtypicai, showing that within the redshift interval displayed, 
each present day MW sized system will have had between 
10 and 100 PopIII halos as progenitors, contributing stars 
with metallicities in —4.0 < [Z/Zq] < —2.4. This means 
that the mass functions of eq. 1 will be sampled densely, 
and the spread between different present day systems will 
be minimal i.e., all MW sized systems should contain the 
same numbers and distribution of stars having metallici- 
ties in the above range. The two curves of Fig. 2 inter- 
sect at a redshift of ~ 15, corresponding to a metallicity 
of [Z/Zq] = —4.3 in the Gnedin & Ostriker picture. At this 



limit many present day MW sized systems will have had 
only one PopIII halo coming from beyond this redshift, and 
many will have had none at all. In this way we can expect 
that the numbers of stars in the Milky Way with metallic- 
ities below [Z/Zq] = —4.3 should be minimal, indeed, no 
such stars have been found to date. 

For even lower metallicities, one goes to larger limiting 
redshifts to reach z ~ 30, Z/Zq ~ and Mt yp i ca i ~ 10 6 Mq, 
the objects responsible for the appearance of the very first 
stars and of starting the reionization of the universe. How- 
ever, as commented above, the chance of finding any such 
object as part of the MW is slim, they will be found as part 
of larger systems, perhaps as intracluster stars. For simplic- 
ity though, we will use the term first stars to identify the 
earliest stars of our Galaxy. 



3 FORMING THE FIRST STARS 
3.1 Disk Formation 

We model each of the PopIII halos discussed above as a small 
galactic system. Systems with masses larger than M cr u cool 
and the gas forms a disk in centrifugal equilibrium (often 
seen in the numerical simulations quoted above), with an 
exponential surface density profile fixed by the total bary- 
onic mass and angular momentum. This last is given by the 
A parameter, which is chosen at random from the distribu- 
tion: 



P(X) = 



<x A (2 T )i/2 



exp 



ln 2 (A/ < A >) 
2al 



T' 



(3) 



with < A >= 0.05 and a\ = 1.0, e.g. Dalcanton et al. (1997) 
and references therein. The initial dark halo profile is cho- 
sen such to have a central constant density core, as seen 
in present day dark halos from dwarf to cluster scales, Fir- 
mani et al. (2000). The formation of the disk modifies the 
gravitational potential in the central regions, we include the 
reaction on the dark halo to this process through an adia- 
batic invariance hypothesis (Flores et al. 1993; Hernandez 
& Gilmore 1998). 

Fig. 3 shows a typical PopIII proto-galaxy, having a 
total mass of 3 x 10 s M@, with a close-to-average angular 
momentum, A = 0.035. The final rotation curve is given 
by the solid curve, with the dotted line giving the surface 
density profile of the disk. 



3.2 Disk Instabilities 

The final temperature of the gas in the disk is taken to 
be 300K, as appropriate for gas where the main cooling 
mechanism is H2 (Ferrara 1998, Susa & Umemura 2000, 
Machacek et al. 2000). This temperature determines c s , the 
sound speed within the disk, which together with re, the 
epicycle frequency obtained from the detailed final rotation 
curve, yields the Toomre's stability parameter for the disk: 



Q 



(4) 



In the regions where Q > 1 the total tidal shears, to- 
gether with the velocity dispersion, are sufficient to stabilize 
the disk locally against its self gravity. These regions are 
therefore not subject to star formation. The dashed curve 
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Figure 3. A typical PopIII galaxy, total mass 3 X 10 Mq, with a 
close to average angular momentum, A = 0.035. The solid curve 
gives the rotation curve. The red line represents the baryonic disk 
surface density profile; together with the estimated sound speed 
in the disk, it gives the inverse of Toomre parameter, Q _1 , shown 
by the dashed curve. The Q = 1 instability threshold is shown by 
the dotted line, with the intersection of this and the Q~ 1 curve 
defining the stable and unstable regions. 




Figure 4. Same as Fig. 3, but for a higher value of the angular 
momentum, A = 0.08. The changes in surface density profile, 
rotation curve (through the gravitational coupling of the disk and 
halo) and Toomre's parameter result in a larger instability radius, 
but a much reduced fraction being susceptible to star formation. 



in Fig. 3 gives the profile of Q~ l in the disk, showing the 
regions interior to the intersection with the Q = 1 line to 
be unstable towards star formation, in this case 85% of the 
disk mass. Fig. 4 is analogous to Fig. 3, but for a system 
having a larger angular momentum, A = 0.08. The changes 
in disk surface density and rotation curve imply than now 
only a smaller fraction (60%) of the disk will form stars. A 
value for the gas temperature lower than the 300K we are 
assuming, due to the small presence of heavy elements in 
these objects, would only increase slightly the value of the 
unstable fraction. Given that the unstable fractions for sys- 
tems with typical values of A are already close to unity, this 
detail does not affect our conclusions. 

Finally, a gas to star conversion efficiency /* = 0.05 
is used to turn the total mass within the instability region 
into a final mass of stars, for each object. Although we have 
calibrated this parameter from local observations (Ciardi et 
al. 2000), since there is no firm theoretical prediction for 
any evolution with redshift of it, we have not explored such 
possibility. 

3.3 Cloud Collisions and Collapse 

Although the above instability criterion might be satisfied, 
it does not necessarily follow that the considered PopIII ob- 
jects will be able to form stars. In fact, it is possible that 
cloud-cloud collisions might disrupt the gas clouds formed 
by the gravitational instability before they have time to col- 
lapse (a point originally made by Kashlinsky & Rees 1983). 
To explore this possibility we calculate the time for colli- 
sions between clouds, t co u and for cloud collapse, tff. This 
last will be given by: 

^(A) 1/2 ' (5) 

where p is the density of the gas clouds. The collisional time 
scale is 



(6) 



with M the number density of clouds, and Aj the Jeans 
length in the disk (in general a function of radius) which 
we take as a representative size for the clouds. Taking M = 
f c p/Mj, where f c is the volume fraction of the disk filled 
with clouds, a filling factor, and the Jeans mass 



Mj 



(¥) 



Gp 



3/2 



since the Jeans length is: 
2\ i/a 



A, 



7TC C 

~~Gp 



(7) 



(8) 



(9) 



we can calculate the ratio of the collision time to the collapse 
time by substituting the values above into eqs. 5 and 6 to 
yield, 

tcoii _ 2^9 
*// fc 

Eq. 9 shows that even if the filling factor for the fragmenting 
disk were unity, which is unlikely, the collapse time of the 
clouds is always shorter than the collision time, at all radii, 
and independent of the value of c s . This is a result of having 
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taken a constant value of c s throughout the disk, which is 
a consequence of H2 being the dominant cooling mechanism 
operating in the primordial metallicity gas. As eq. 9 shows, 
for all the scales of our problem, gravitational cloud collapse 
can proceed unimpeded by cloud-cloud collisions, which we 
shall not consider further. 

3.4 The First Stars Today 

Once the total gas mass turned into stars as a function of the 
metallicity has been determined, we use the IMF of Larson 
(1998) to turn this into a total number of stars. This function 
offers a convenient parameterization of the IMF: 

dN/dlogm oc (1 + m/m s ) _1 - 35 . (10) 

In the above equation m s is a characteristic mass scale, of 
order O.35M0 for a present day solar neighbourhood IMF. 
This scale mass can be increased to explore the consequences 
of a top heavy IMF applying to stars of metallicities much 
lower than solar, as the ones we are treating. 

We combine this sampling of the IMF using m s = 
O.35M0 with the most recent Padova stellar models of pri- 
mordial metallicities of Girardi et al. (2000) to exclude those 
stars which due to their age have by now evolved off the tip 
of the RGB, and are hence no longer visible. In this way 
we obtain a sample of present day stars. Sampling the IMF 
from a fixed lower luminosity up to the tip of the RGB of 
the theoretical isochrones allows us to obtain the number of 
stars visible today upwards of a certain limiting luminosity, 
as a function of metallicity. Our final samplings include only 
stars between O.GMq and the tip of the RGB of the Padova 
isochrones (of around I.OMq), which is the range over which 
stellar models were calculated. 

As the total mass in PopIII objects is high, both the 
PopIII halos mass functions and the A distributions are sam- 
pled densely, yielding a low spread for the predicted final 
number of stars. 

To summarize, each of the PopIII object selected from 
the mass function of eq. 1 is carefully modeled as a mini 
galaxy, with a value of A selected at random from eq. 3. The 
resulting system is then treated as a mini galaxy where the 
cooling of the baryons leads to disk formation; subsequently 
a disk instability criterion is used to obtain the total mass 
in stars contributed by each object. Finally, through the 
use of an IMF and new low metallicity stellar models, this 
total mass is turned into a number of stars visible today. 
This is repeated until the total PopIII object mass equals 
that within the bounded region described in Section 2, for a 
given limiting redshift, which corresponds to a fixed maxi- 
mum metallicity. This process is repeated for different limit- 
ing redshifts, to construct a total mass distribution function 
for the low metallicity stars for the present day MW. The 
curve in Fig. 5 shows the number of MW stars with mass in 
the above range and having metallicities lower than a given 
value. 



4 COMPARISON WITH DATA 

To compare our results with observations we take the sample 
of Ryan & Norris (1991) of low metallicity halo stars in 
the local disk. They include a magnitude limited sample, 
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Figure 5. Average number of PopIII stars with metallicity lower 
than [Z/Zq] expected in MW sized galaxies today, for a solar 
neighbourhood IMF (solid curve). The dots with error bars show 
the observations of Ryan & Norris (1991) and Beers et al. (1998), 
assuming an isothermal distribution for the stars in the halo (see 
text). The labeled arrows show the values of m a required to rec- 
oncile model and data, at z = 6.5 and z = 9.0. 

(their NLTT sample) which is complete down to O.GMq and 
V — 13.0, the same limit we use to sample our IMF. From 
the theoretical isochrones, this corresponds to a distance 
of 290 pc, which we take as the radius of the sample. In 
that sample 5 stars were found with [Z/Zq] < —3.0. To 
obtain the total number of stars in the Galactic halo at this 
metallicity, we must introduce an assumption regarding the 
distribution of these stars. 

If the PopIII stellar systems are disrupted while they 
are incorporated into the growing MW, their stars will to- 
day be found almost exclusively in the Galactic halo, while 
if the initial stellar systems form tight clusters, they might 
survive disruption and through dynamical friction end up in 
the galactic bulge. Which of these two processes dominates 
will depend on the mass ratio of the mergers through which 
the PopIII systems reach the MW. If these mass ratios are 
close to one, disruption will be important, while if they are 
accreted late into a fully formed MW PopIII objects might 
remain bound. From N-body simulations White & Springel 
(2000) have argued that most of the lowest metallicity stars 
in a present day galaxy will be found in the stellar halo, 
and not in the bulge. We therefore make the assumption 
that the lowest metallicity stars share the isothermal dis- 
tribution of the dark halo, out to a distance of 100 kpc, 
representative of the total extent of the halo. The actual 
density profile of these objects might more closely resemble 
that of the old globular cluster population, having an R~ 3 
profile. Fortunately, fixing the normalization point at the 
solar radius (8.5 kpc) yields total numbers which are largely 
independent of this change in the assumed profile, e.g. for 
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Figure 6. Inferred values of m s as a function of redshift (points 
with error bars). The solid curve gives the redshift evolution of 
the Jeans mass of cold star forming clouds, identified with m s , 
resulting from the temperature evolution of the cosmic microwave 
background (CMB). 



Figure 7. HR diagram for 2286 extremely low metallicity stars. 
Accurate isochrones for these stars are used to obtain good es- 
timates of their present observational characteristics, which can 
be used to provide guidelines of their current photometries and 
hence aid detection strategies. 



core radii in the R~ 3 profile going from 0.2 to 2 kpc (a core 
radius must be introduced in the R~ 3 profile to avoid the 
central divergence), the total numbers of stars change by 
a factor 1.24 to 0.34, respectively, small compared to our 
error bars in figures 5 and 6. It is worth pointing out that 
the assumption of a particular density profile for the current 
low metallicity halo stars is not intrinsic to our results, but 
becomes necessary only for data comparison purposes. 

Taking the Ryan & Norris (1991) data plus the hypoth- 
esis of an isothermal distribution discussed above, we de- 
rive the total number of stars with [Z/Zq] < —3.0 in the 
Galaxy as 4.43 x 10 6 . To obtain numbers at other metallic- 
ities we use the relative metallicity distribution function of 
the more recent Beers et al. (1998) work, which combines 
several new samples. These observations are shown in Fig. 5 
by the points with error bars, which reflect the uncertainties 
in the normalization. 

The stellar halo might not extend as far out as 100 kpc, 
but some PopIII stars might be found in the bulge, so the 
interpretation of data used in Fig. 5 is probably a lower limit. 
On the other hand, given the uncertainties in the primordial 
IMF, our prediction is a strict upper limit. 

It is clear that our results fall somewhat above the ob- 
servational measurements, implying that our assumption of 
a constant m 3 = O.35M0, essentially a present day solar 
neighbourhood IMF, is not valid. Reconciling our estimates 
with the observations requires not only the assumption of 
a higher m s in the past, but throughout the redshift range 
studied, an increasing trend for this value. Note that our 
derived numbers of PopIII stars scale linearly with /* 

Several authors, e.g. Larson (1998) have pointed out 



that the increase of the CMB temperature with redshift 
alone, implies an increase in the Jeans mass for the star 
forming clouds, once one is beyond z = 2, where the CMB 
temperature equals the « 8 K of cold star forming clouds. 
It has therefore been suggested that at high redshift the 
IMF should be increasingly weighted towards larger stellar 
masses. The labeled arrows in Fig 5. show the values of m 3 
which we require at z = 6.5 and z = 9.0, to match the 
observed points. 

In Fig. 6 we show the values of m s which we infer for the 
IMF of PopIII objects, by imposing that model and observa- 
tions agree. The solid curve indicates the redshift evolution 
of the Jeans mass in star forming clouds, which we identify 
with m s . Although for z < 9 our results are consistent with 
the increase in m 3 coming solely from that in the Jeans mass 
due to the CMB, the trend for 6 < z < 9 is clearly for a rise 
in m s beyond what this effect can account for. In fact, as 
Larson (1998) points out, this is only the minimum increase 
to be expected, as other heating sources, e.g. massive stars 
themselves or a UV background, were likely to be present. 
Furthermore, a reduced opacity consistent with these low 
metal abundances would have the effect of skewing the IMF 
towards larger masses, thus operating in the same direction. 

For z > 9, the values of the characteristic mass appear 
to stabilize around the value m s = 10 — 15Mq , indicative of 
a very top heavy IMF at these early epochs, although the 
increase of our error bars makes determining trends in this 
redshift range harder. 

Finally, Fig. 7 presents an HR diagram containing over 
2000 stars, produced from the low metallicity stellar models 
used, convolved with errors appropriate to HST observations 
of halo stars (Gallart et al. 1999; Hernandez et al. 2000), for 
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m a = 5.0Mq. The inclusion of detailed atmosphere calcula- 
tions and spatial information to yield detailed observational 
predictions will form part of a forthcoming paper. 



5 SUMMARY 

From our study we conclude that: 

• If the standard hierarchical merger scenario is correct, 
in the MW one should not expect to find stars having metal- 
licities lower than [Z/Zq] = —4.3, consistent with current 
observational limits. 

• We infer the IMF of PopIII stars to have been strongly 
weighted towards high masses, and increasingly more so in 
the range 6 < z < 9, this levelling off at z ^ 9 at a charac- 
teristic mass scale m a = 10 — 15Mq. 

• Observational limits on the numbers and metallicity 
distribution of very low metallicity stars in the MW can 
place sharp restrictions on the IMF of the first stars. This 
can be achieved through detailed modelling of the full range 
of phenomena giving rise to such stars, of the type attempted 
here, thus necessarily establishing a connection with the un- 
derlying structure formation scenario. 
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